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EXPANSIONS  OF  CUMULATIVE  DISTRIBUTION  FUNCTIONS 
DIRECTLY  FROM  CHARACTERISTIC  FUNCTIONS 

INTRODUCTION 

The  performance  analysis  of  nonlinear  signal  processors  in 
noise  can  often  be  accomplished  in  closed  form  in  terms  of  the 
characteristic  function  of  the  system  output  decision  variable; 
for  example,  see  [1  -  4].  However,  due  to  analytical  complexity, 
determination  of  the  corresponding  probability  density  function 
and  the  cumulative  or  exceedance  distribution  functions  then 
generally  requires  some  sort  of  numerical  procedure,  such  as  a 
fast  Fourier  transform;  for  example,  see  [5], 

Alternatively,  expansions  of  the  probability  density  function 
or  cumulative  distribution  function  in  a  Hermite  or  generalized 
Laguerre  series  is  possible  if  the  system  output  high-order 
moments  or  cumulants  can  be  easily  determined  [1,6],  However, 
this  approach  of  breaking  a  given  characteristic  function  down 
into  its  moments  or  cumulants,  followed  by  a  build-up  of  a  series 
expansion  from  coefficients  determined  recursively  from  the 
moments  or  cumulants  [6],  can  involve  a  large  amount  of  number¬ 
crunching  with  its  attendant  round-off  error.  What  is  desirable, 
instead,  is  a  method  of  proceeding  directly  from  the  closed  form 
characteristic  function  to  series  expansions  for  the  probability 
density  function  and  cumulative  distribution  function.  A  step  in 
this  direction  was  taken  recently  in  [3,4],  where  the  receiver 
output  characteristic  function,  for  a  multiple-pulse  processor  in 
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a  noisy  medium  with  correlated-fading,  was  expanded  in  a  series 
involving  only  positive  terms  and  coefficients.  Here,  we  will 
generalize  those  results  and  obtain  an  extension  of  the 
generalized  Laguerre  expansion. 

The  ideal  goals  of  this  analysis  would  be  as  follows: 
the  new  series  expansions  should  be  rapidly  convergent  and 
involve  as  few  terms  as  possible;  the  series  should  contain  some 
adjustable  parameters  that  can  be  modified  to  achieve  more  rapid 
convergence,  especially  for  difficult  situations  such  as  very  low 
false-alarm  probabilities  or  very  high  signal-to-noise  ratios; 
the  individual  terms  in  the  series  should  be  easily  computed, 
preferably  by  efficient  recursions;  the  series  for  the 
probability  density  function  should  be  easily  integrable  to  yield 
the  cumulative  or  exceedance  distribution  functions;  and  the 
series  should  involve  only  positive  terms,  if  feasible,  in  order 
to  retain  high  accuracy  in  the  tail  probability  regions. 

The  following  results  have  achieved  partial  success  in 
realizing  these  goals.  The  expansions  are  not  limited  to 
orthogonal  expansions,  such  as  the  Laguerre  and  Hermite  cases. 
There  exists  the  possibility  of  creating  additional  convergence 
factors  in  the  expansions  which  speed  convergence  significantly. 
And,  arbitrary  signal-to-noise  ratios  and  thresholds  can  be 
easily  accommodated  without  sacrificing  accuracy.  There  still 
exists  the  problem  of  deciding  how  best  to  pick  the  free 
parameters  in  the  expansions  for  most  rapid  convergence.  It  is 
generally  not  simply  to  make  a  few  low-order  coefficients  zero; 
this  is  consistent  with  [6;  page  3  and  examples]. 
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EXPANSION  OF  CHARACTERISTIC  FUNCTION 


Before  presenting  the  new  expansions,  it  is  worthwhile  to 
review  some  earlier  results.  The  generalized  Laguerre  expansions 
of  a  probability  density  function  p(u),  for  a  positive  random 
variable,  and  its  corresponding  characteristic  function  (Fourier 
transform)  f(K)  are  given,  respectively,  by  [6;  (102)]  as 

P(u)  =  exp(-u/P)  ^  L(v-l)(u]  foru>0,  (1) 

pv  T(\>)  n=0  n  n 


fU) 


1 

(1  -  UP)V 


00 


n=0 


(\>) 


n 


n  n ! 


(2) 


Except  for  the  leading  factor  1/(1  -  i£P)v,  the  basic  idea  here 
can  be  interpreted  as  the  expansion  of  the  characteristic 
function  f(l)  in  a  power  series  of  the  quantity  -if;/(l  -  i£P). 
There  are  two  free  parameters  in  (2),  namely  v  and  p,  which  can 
be  chosen  to  speed  up  convergence,  that  is,  hasten  the  decay  of 
the  terms  in  (1)  and  (2).  Observe  that  an  =  an(v,P). 

An  alternative  approach  is  given  by  the  following  expansion 
of  the  characteristic  function,  involving  only  powers  of  the 
factor  1/(1  -  i£p);  see  [3;  appendix  D]  and  [4;  appendix  A]: 
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Again,  there  are  two  free  parameters,  v  and  p,  and  an  =  a  (v,P). 
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LINEAR  FRACTIONAL  TRANSFORMATION 

Both  of  the  expansions  in  (2)  and  (4),  above,  are  special 
cases  of  a  linear  fractional  (or  bilinear)  transformation 

z  -  f  I  ifg  -  iE  -  7^-ik  '  Y  *  pa  ,  e>0,  (5) 

between  the  complex  £  and  z  planes.  A  linear  fractional 
transformation  is  the  most  general  analytic  transformation 
yielding  a  one-to-one  correspondence  between  all  points  in  the 
simple  l  and  z  planes  [7;  page  365].  (The  transformation  in 
(2)  corresponds  to  a  =  0,  y  =  p,  while  that  in  (4)  utilizes 
y  =  0,  a  =  1.  ) 

We  observe  the  following  associations: 

£  =  0  maps  into  z  =  a  ; 

£  =  e  <<  1  maps  into  z  ~  (a  -  iey)(l  +  iE0)  ~  a  +  ie(0a  -  y)  ; 

£  =  °°  maps  into  z  =  y/0  •  (6) 

Then,  for  positive  parameters  a,0,y,  the  upper  half  S-plane  maps 
into  the  interior  of  a  circle  in  the  z-plane,  with  center  c  and 
radius  r',  where 


This  result  holds  regardless  of  whether  0a  >  y  or  (3a  <  y. 
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EXPANSION  PROCEDURE 

To  demonstrate  how  the  expansions  are  developed,  suppose 
characteristic  function  f(l)  is  given,  and  that  its  asymptotic 
behavior  is  governed  by 

f ( £ )  ~  c/(U)v  as  l  -»  «  ,  v  >  0  .  (8) 

This  behavior  determines  parameter  \> .  Now,  define  auxiliary 
function 

*<z>  *  f(£  -17^-fe)  •  <9> 

Form  the  following  product,  and  then  develop  the  product  in  a 
power  series  expansion  in  z  according  to 


n 


n 


n 


=  fn(v,p,a,y) 


(10) 


Then,  from  (9)  and  (5),  the  original  characteristic  function  is 


f(0 


-  a  ~ 

1  -  UpJ 


1 

(1  -  HP)V 


(11) 


The  scale  factor  y  -  is  used  in  the  product  on  the  left  side 
of  (10)  in  order  to  get  the  more  compact  form  (11).  The 
asymptotic  behavior  of  expansion  (11)  agrees  with  (8),  due  to  the 
introduction  of  leading  factor  (y  -  0a)/(y  -  0z)  =  1  -  i£0  in 
(10),  before  the  power  series  expansion  in  z  is  accomplished. 

The  coefficients  { f n |  in  (10)  and  (11)  are  functions  of  the 
four  parameters  v,p,a,y.  If  a  =  0,  the  quantity  fn(v,(3,0,y)  yn 


5 


TR  10371 


in  (11)  can  be  denoted  as  fn(v,(3),  giving  form  (2).  On  the  other 
hand,  if  y  =  0,  the  quantity  f  (v,p,a,0)  an  in  (11)  can  be 
denoted  as  fn(v,0),  giving  form  (4). 

Although  four  parameters  are  available  in  (11),  only  three 
are  fundamental;  that  may  be  seen  by  writing  (11)  in  the  form 


f(S) 
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(1  -  UP) 
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00 


n=0 
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n  -  uv/ar1 

1  -  Up  ' 


(12) 


provided  that  a  ^  0.  Coefficient  f  an  is  then  a  function  only 
of  the  three  parameters  v,  ft,  a/y.  Since  characteristic  function 
value  f(0)  =  1,  there  follows  from  (12)  the  useful  constraint 


00 


n=0 


n 


1  . 


(13) 


Another  possible  triplet  of  parameters  is 

v  asymptotic  decay  ; 

0  scaling  constant  ; 

(3a/y  matching  constant  .  (14) 

Whichever  set  of  three  parameters  are  used,  they  should  be  chosen 
so  that  series  (11)  or  (12)  has  rapidly  decaying  coefficients 
I f n } ,  if  possible. 
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EXPANSION  OF  PROBABILITY  DENSITY  FUNCTION 


In  this  section,  the  probability  density  function 
corresponding  to  characteristic  function  expansion  (11) 
will  be  developed.  For  future  use,  we  define 


r  -  ,  V  *  0  . 


(15) 


We  begin  by  observing  that  the  numerator  term  (a  -  Uy)n  in 
(11)  can  be  expressed  as 


<«  -  «v)n  - 


+  1  -  UP 


-  &)"£(:)  ^  <i  -  iwB- • 


(16) 


This  enables  us  to  express  the  n-th  term  in  (11)  according  to 

Trris? (H-Hfl" -  SI" SO'- » 


Then,  characteristic  function  (11)  takes  the  desired  form 


£<E>-Sfn 


-v-m 


(18) 


The  expansion  of  probability  density  function  p(u), 
corresponding  to  this  form  of  f(l)  in  (18),  is 


p(u)  = 


;u/p)v~1  exp(-u/ 
P  T(v) 


SnH'")  ^r  u>  0,(19) 
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where  auxiliary  function 


Sn(*,v) 


jFj(-n;v;-x| 


n 


<v>n 


(20) 


Here,  we  used  [8;  (13.6.9)].  The  recursion  for  Sn(x,v)  is,  by 
means  of  [8;  (13.4.1)], 


Sn(x,v) 


( v+2n-2+x )  Sn_1(x,v)  -  (n-1)  Sn_2(x,v) 

v+n-1 


(21) 


for  nil,  with  Sq(x,v)  =  1.  Value  S_^(x,v)  can  be  arbitrarily 
set  to  zero,  since  it  gets  multiplied  by  zero  when  n  =  1  in  (21). 

From  (15),  if  y  <  then  r  >  0  and  all  the  terms  in 

Sn(ru/0,v)  in  (20)  are  positive.  Sequence  jSn(ru/p,v)}  in  (19) 
does  not  decay  rapidly  to  zero  with  n;  in  fact,  the  sequence  can 
grow  slowly  with  n,  depending  on  the  polarity  of  r.  Convergence 
of  (19)  must  be  furnished  by  sequence  fn  (y/p)n. 

We  could  express  the  probability  density  function  in  (19)  as 


P(u) 


(u/P)v  1  exp(-u/P) 
P  F(v) 


00 


n=0 


n 


(v-1) f^ru) 
n  <5  J 


(22) 


(Then,  by  using  [9;  7.414  8],  the  corresponding  characteristic 
function  is  readily  obtained  as  (11),  which  constitutes  a  check 
on  the  above  developments.)  Notice  that  the  argument,  -ru/p, 
of  the  Laguerre  polynomial  Ln  in  (22)  is  not  the  same  as  the 
argument,  u/p,  of  the  exponential,  unless  r  =  -1,  that  is, 
unless  we  are  dealing  with  the  special  case  of  a  =  0. 
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For  a  =  0,  the  results  above  reduce  to 
,  v— 1 


P(U)  =  ZZ  fn  (^)n  ^  ^(f)  for  u>0  (23) 


t(l)  = 


(1  -  U(3)v  n=0 


op  /  ■  f  \  n 

=  £n(i^fe)  ' 


(24 


These  agree  with  the  usual  generalized  Laguerre  expansion  in 
[6;  (102)]  when  we  set  y  =  |3  and  identify  v,  f  here  with 
a+1,  an  (a+l)n/n!  there,  respectively. 

On  the  other  hand,  in  the  limit  as  y  •*  0,  observe  from 
(20)  that 

»)“ -.W  ■  ft)“  S 13  ik  l1^  t)“  ■ 
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1  f* 

xu)n 

IpJ  (>»)n 

l  Y 

pJ 

(v)n  l 

Pi 

(25) 


Then,  the  result  in  (19)  reduces  (as  y  ->  0)  to 

v— 1 


PCI  -  (U/P17-r^?i'U^1  EZ  f.  for  u  >  0  ,  (26| 


n=0  n  (v>n 


while  (11)  becomes 


£<£l  = 


n 


(27) 
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These  agree  with  [3;  ( D— 1 5 )  and  (D-13)],  respectively,  when  we 
set  a  =  1,  £  =  2,  and  identify  v,  f  here  with  K,  F  gn  there, 
respectively. 

Relations  (22)  and  (19)  are  extensions  of  the  generalized 
Laguerre  expansion  of  a  probability  density  function.  These  are 
the  desired  results,  in  so  far  as  the  probability  density 
function  is  concerned;  we  now  proceed  to  expansions  for  the 
corresponding  cumulative  distribution  function. 
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EXPANSION  OF  CUMULATIVE  DISTRIBUTION  FUNCTION 


The  cumulative  distribution  function  is 

v 

C(v)  =  J  du  p(u)  .  (28) 

—  00 

We  will  obtain  expansions  for  C(v)  in  two  different  ways,  each 
with  its  own  merits  and  limitations. 


FOURIER  TRANSFORM  OF  C(v) 

The  Fourier  transform  of  the  cumulative  distribution  function 
C(v)  is  considered  in  appendix  A;  it  is  given  in  terms  of  the 
characteristic  function  f(£)  by 

g(£)  *  |  dv  exp{ i£v )  C(v)  -  ^  f(l)  for  0  <  <  a  ,  (29) 

where  is  the  imaginary  part  of  £.  Constant  a  is  positive,  and 

can  be  +°°,  as  for  a  random  variable  which  is  always  positive. 
Then,  from  (8),  the  asymptotic  behavior  of  g(S)  is 


g(U  ~ 


-c 


<U) 


\>+l 


as  l  ->  00  . 


(30) 


Analogous  to  (9)  and  sequel,  we  define  auxiliary  function 


a<z>  5  9^  -  T  ■ 


_  Y  -  Pz 
a  -  z 


f(z)  / 


(31) 


where  we  used  (29).  Now,  we  form  the  following  product,  G(z), 
and  develop  the  product  in  a  power  series  in  z  according  to 
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(  -  R  \  V+*  00 

;(2>  s  2(z|  =  £Z  9n  z"  ;  9n  =  9n(v,p,o, 


V)  •  (32) 


In  order  to  find  coefficients  (gn),  we  use  (31)  and  (10)  to  find 


G( z )  =  - 


(y  -  pz)  (a  -  z) 


f(z)  =  C  f  z“  .  (33) 

— '  '  a  -  z  1 — *  n  v  ' 

n=0 


Now,  if  a  /  0,  we  expand  (a  -  z)-*  =  a-*  (1  -  z/a)~*  in  a 
power  series  in  z,  and  combine  the  two  series  in  (33),  thereby 
obtaining  coefficients  (gn)  in  (32)  as 


gn  -  El  f„  om  -  i  9n  ,  +  *„  for  n  2  0  ,  (34) 

Jn  n+l  * — s*  m  a  ^n-l  an 

a  m=0 


with  g =  0.  Thus,  coefficients  ign)  are  readily  found  in  terms 
of  the  earlier  coefficients  { f n ]  in  (10),  at  least  when  a  ^  0. 
Finally,  using  (31)  and  (32),  we  have  the  desired  expansion  for 
the  Fourier  transform  (29),  namely 

9<f>  -  a(z  “  T-TT^)  =  ‘  v+1  ^  9n  (H-lfp)  '  <35> 


where  gn  =  gn(v,p,a,Y). 
Observe  from  (34)  that 


lim  g  an+1  =  ( pa  -  y)  2IZ  f  am  =  pa  -  y  , 

n-»°>  n  m=0  m 


(36) 


where  we  used  (13).  Therefore,  the  asymptotic  decay  of 
coefficients  jgn)  is  according  to 
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9 


n 


-  V 
an+1 


as  n  ->  «  . 


(37) 


This  result  demonstrates  that  series  (35)  will  converge  if  |y|  is 
less  than  |£a|. 

Then,  by  reference  to  the  Fourier  transform  pair  (11)  and 
(19),  the  cumulative  distribution  function  (inverse  Fourier 
transform)  corresponding  to  (35)  is  immediately  available  as 


C(  v) 


=  (v/P)  exp(-v/g) 
13  r(v+i) 


n=0 


(fl‘  »nH^)  - 


where  we  also  used  (20).  These  results  apply  only  for  a  ^  0. 

The  series  converges  if  |y|  <  |0a|,  as  may  be  seen  by  reference 
to  (37)  and  [1°/  (8.22.1)].  The  relations  in  (38)  are  the  main 
results  of  this  subsection. 

If  0  <  y  <  Pet/  then  r  >  0  from  (15),  and  all  the  terms  in 
Sn( rv/p, v+l )  are  positive.  Again,  notice  that  the  argument, 
-rv/p,  of  the  Laguerre  polynomial  is  not  the  same  as  the 
argument,  v/p,  of  the  exponential,  unless  r  =  -1,  that  is,  unless 
we  are  dealing  with  the  special  case  of  a  =  0.  However,  this 
special  case  cannot  be  obtained  from  the  above  by  simply  letting 
a  -»  0  in  (38);  see  the  comment  under  (33),  where  property  a  ^  0 
was  explicitly  used  in  developing  the  series  expansion  for 
(1  -  z/a)'1. 
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SPECIAL  CASE  OF  a  =  0 


From  (32)  and  (33),  when  a  =  0,  we  have 


2(z)  =  [- 


v+1  » 


G(z)  = 


*  n- 1 

fn  z 


(39) 


with  f  =  fn(v,0,O,Y)*  Then,  (35)  yields 


>’  ■  *[■  ■  h%1  ■  ;rrS7rT  5  <■  ‘  • 


Again,  we  could  let  fn(v,0,O,y)  y  =  fn(v,P),  where  the  latter 
quantity  is  independent  of  y. 

The  cumulative  distribution  function  term  corresponding  to 
the  n  =  0  term  in  (40)  requires  special  treatment;  we  have 


q  (£\  =  - 1 -  f  i — 

°  (1  -  i^)V+1  ° 


-U(l  -  UP) 


(41) 


Since  we  know  the  Fourier  transform  correspondence 


(1  -  UP)' 


v-1 

u  exp(-t 
PV  r(v) 


for  u  >  0  , 


(42) 


the  9q(1)  term  corresponds  to  the  cumulative  distribution 
function  component  (see  (28)  and  (29)) 


co<v>  “  £o  I 


V—  1  / 

u  exp(-u/( 
T(v) 


■o  T(v) 
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-  £0  <V/P)^,;;p<-V/P>  iri(„ltv;v)  for  v  >  0  ,  ,43, 

where  y(v,x)  is  the  incomplete  gamma  function  [8;  (6.5.2)],  and 
jFj  is  the  confluent  hypergeometric  function  [8;  (6.5.12)]. 

The  remainder,  g^(S),  of  g(£)  in  (40)  can  be  inverse 
transformed  immediately  by  comparing  Fourier  transform  pair  (11) 
and  (19).  Its  cumulative  distribution  function  component  is 


Cx(v) 


■’nM' 


(44) 


for  v  >  0.  The  total  cumulative  distribution  function  is  given 
by  the  sum  of  (43)  and  (44): 


C( v)  =  CQ(v)  +  Cj ( v)  = 


(45) 


r(v+l) 


[ f o  lFl(1;1+v?p)  +  fn  (p)  Sn-1  (  P'v+1)l  ' 


This  expansion  for  the  cumulative  distribution  function  applies 
for  a  =  0.  Again,  we  could  let  fn(v,p,0,y)  yn  =  ?n(v,p). 

An  alternative  expression  for  C(v)  is  available  when  we  use 
identity  (20)  on  result  (45).  When  this  is  done,  the  result  for 
C(v)  in  (45)  agrees  with  [6;  (95)],  when  we  set  y  =  (3  and 
identify  v,  fn  here  with  a+1,  an  (a+l)n/n!  there,  respectively. 
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SPECIAL  CASE  OF  y  =  0 

As  another  check  on  general  result  (38),  we  consider  letting 
Y  -»  0.  From  (25),  we  have 

*)"  Mir'H  *  T^n;  H)n  «■  r  *  o .  («> 

Then,  from  (38),  when  y  =  0,  the  cumulative  distribution  function 
becomes 


C(  v) 


(v/p)v  exp( -  v  y““’  a  igy/gl-  = 

p  r(M+l)  yn  (v+l)n 


=  (v/p)  exp( -v/p ) 
r( v+i > 


00  ,  n  » 

EZ  EZ  £m  “m 

n=0  vm=0  ra  ‘ 


(y/fi) 

(v+l) 


n 


n 


for  v  >  0 


(47) 


rri 

Here,  also,  we  could  let  fm(v,p,a,0)  a  =  fm(v,P). 

The  quantity  1  -  C(v)  is  a  rearrangement  of  [3;  ( D—  21)]; 
agreement  can  be  confirmed  when  we  set  a  =  1,  p  =  2,  use  (34), 
and  identify  v,  fn  here  with  K,  F  gn  there,  respectively. 
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INTEGRATION  OF  PROBABILITY  DENSITY  FUNCTION  EXPANSION 

At  the  beginning  of  this  section,  it  was  noted  that 
expansions  for  the  cumulative  distribution  function  C(v)  would  be 
obtained  by  two  different  methods.  The  results  above  were 
obtained  by  means  of  the  Fourier  transform  of  C(v);  now,  we 
integrate  directly  on  the  probability  density  function  p(u). 

From  (12),  the  characteristic  function  f(£) ,  for  arbitrary 
v,  0,  a,  Y/  can  be  expressed  as 

00 

f(U  =  YU  fn  a"  F  U)  ,  a  *  0  ,  (48) 

n=0  n  n 


where  characteristic  function  component 


Fn(E) 


(1  -  Uv/cOn 
(1  -  UP)v+n  ' 


v  >  0  . 


(49) 


Component  ?n(Z)  is  a  characteristic  function,  in  that  Fn(0)  =  1 
for  all  n.  Then,  from  (19),  the  probability  density  function 
corresponding  to  (48)  is 

00 

P(u)  =  ZZ  fn  Qn  pn(u>  for  U  >  0  ,  (50) 

n=0  n  n 


where  probability  density  component 


pn(u) 


f_^|n  (u/0)v~1  exp(-u/p) 
UaJ  0  r(v) 


.  fru 

’ni  0,v 


for  u  >  0  . 


(51) 


Component  pn(u)  is  a  probability  density  function,  in  that 
S  du  pn(u)  =  1  for  all  n. 
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see  appendix  B  for  additional  useful  recurrence  relations  on  this 
incomplete  gamma  function.  Component  Cn(v)  is  a  cumulative 
distribution  function,  in  that  Cn(°°)  =  1  (using  C(a,°°)  =  1). 
Relations  (52)  -  (53)  are  the  main  results  of  this  subsection. 

The  cumulants  corresponding  to  an  individual  component  may  be 
found  readily  from  (49);  the  k-th  order  cumulant  is 

Xn(k)  =  ( k-1 )  I  ((v+n)  Pk  -  n  (*)  j  for  k  M  .  (55) 

In  particular, 

Xn(l)  =  (v+n)  P  -  n  J  , 
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Xn( 2 )  =  (v+n)  p2  -  n  ^  . 
n  a* 

Therefore,  if  y  <  pa, 


and 


Xn(  1 )  -*  00  as  n  -»  «  , 
Xn(2)h 

^nr  * 0  as  n *  - 


(56) 


(57) 


(58) 


That  is,  the  mean  of  probability  density  function  component  Pn(u) 
tends  to  °»,  while  the  component  concentrates  around  its  mean. 

This  implies  that,  for  fixed  threshold  v,  the  cumulative 
distribution  function  component 

Cn(v)  0  as  n  •*  ®  ,  (59) 

provided  that  y  <  pa.  Thus,  expansion  (52)  has  two  convergence 
factors  when  y  <  pa,  but  otherwise  arbitrary  p,  a,  y;  recall 
(13).  Also,  if  y  <  pa,  then  r  >  0,  and  all  the  terms  in 
component  Cn(v)  in  (53)  are  positive. 
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SPECIAL  CASE  OF  a  =  0 


Instead  of  letting  a  ->  0  in  (52)  and  (53),  we  return  to  (23), 
which  gives  the  probability  density  function  p(u)  for  a  =  0.  The 
corresponding  cumulative  distribution  function  is 


for  v  >  0,  where  we  used  [6;  appendix  C).  The  result  in  (60) 
agrees  with  [6;  (91)  and  (95) J.  If  desired,  we  could  again  let 
f n ( v ' ^ , 0 ' Y )  Yn  =  fn(^/P)/  indicating  a  two-parameter  dependence. 

SPECIAL  CASE  OF  y  =  0 

Observe,  from  (53),  that 

cn(v)  ■*  c(v+n,^)  as  y  ->  0  ,  (61) 

meaning  that  we  have,  from  (52),  cumulative  distribution  function 

C ( v)  -»  fn  ^  2  v+n,^  as  Y  ■>  0  .  (62) 

n=0  n  1  p' 

The  corresponding  limiting  exceedance  distribution  function, 

E(v)  =  1  -  C(v),  agrees  with  (3;  (D-21)],  when  we  set  a  =  1, 

0  =  2,  and  identify  \>,  fn  here  with  K,  F  gn  there,  respectively, 
and  recognize  that  C(n,x)  here  is  1  -  H  j.(x)  there.  Also,  since 
C(°°,x)  -  0  according  to  (B-4),  the  sum  in  (62)  has  two 
convergence  factors  as  n  Also,  f  (v,0,a,O)  an  =  ?n(v,P). 
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TYPICAL  CHARACTERISTIC  FUNCTION  COMPONENTS 

Up  to  this  point,  the  given  characteristic  function  f(£)  has 
been  general;  the  power  series  expansion  developed  in  (10) 
allowed  for  arbitrary  f(S)  and  f(z)  in  (9).  Therefore,  all  the 
expansions  developed  in  the  previous  three  sections  had  general 
coefficients  { f n }  and  { gn }  in  the  expansions;  see  (38)  and  (52), 
for  example. 

Now,  we  will  consider  some  typical  characteristic  functions 
that  have  occurred  in  recent  signal  processing  studies  [3,4],  and 
will  develop  the  specific  equations  for  the  coefficients  { f n ) . 

The  four  types  of  characteristic  function  components,  that  we 
will  concentrate  on,  are  labeled  as  follows: 


A: 

f(S) 

-  (1  - 

-  Uafv  , 

(63) 

B: 

f(S) 

M 

=  rri 

m=l 

\  -v_ 

l1-^)  - 

(64) 

C: 

fU) 

=  exp 

[h  £  ,  )  . 

m=l  1  lt,cm 

(65) 

D: 

f(U 

=  exp 

r  M  6m  \ 

fe  >  -  '«J  - 

(66) 

These  individual  forms  may  not  be  valid  characteristic  functions 
for  some  choices  of  parameter  values;  that  is,  the  corresponding 
probability  density  functions  may  go  negative  somewhere. 
Nevertheless,  these  functions  frequently  occur  in  combinations 
and  are  useful  as  building  blocks  for  more  complicated  forms. 
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TYPE  A,  (63) 


The  characteristic  function  of  interest  is  given  in  (63). 

Ihe  asymptotic  decay  is  in  agreement  with  (8).  Function  f_(z)  in 
(9)  is  then 


f  (z)  = 


v  -  pz _ r 

y-aa-z((5-a)J 


(67) 


The  produce  function  in  (10), 

F(Z)  =  fy-~4a)V  f (z)  =  YZ2  f  zn 

v  '  ly  -  0zj  ^  n 

takes  the  form 


F ( z )  =  t  (1  -  hz)  , 


(69) 


for  y  ^  aa. 

Probably,  the  most  efficient  way  to  determine  the 
coefficients  { f n |  in  expansion  (68)  is  to  take  a  logarithm 
of  (69)  and  expand  in  a  power  series: 


In  F ( z )  =  In  t  -  v  ln(l  -  hz)  =  )  \  e  zn  ,  (70) 

n=0 


where  the  latter  coefficients  are  given  by 


e  =  In  t  ,  e  =  —  hn  for  nil. 
o  n  n 


(71) 


Then,  exponentiation  of  (70)  immediately  leads  to  coefficients 


n  , 

f  =  t  ,  f  =  -  hK  f  .  for  n  i  1  . 
o  '  n  n  n-k 


(72) 
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Here,  we  used  the  following  useful  recursion  [6;  appendix  A]: 


00  00 

’(z)  =  exp(]T^  e  zn]  =  []  f  zn 

11  /  n=n  11 


according  to  (70)  and  (68),  with 


f  =  exp(e_ )  ,  fn  =  -  >  k  e,  f  ,  for  nil. 
o  o  n  n  k  n-K 


Result  (72)  is  an  efficient  recursion  for  the  desired  coeffi¬ 
cients  { f n } .  If  p  =  a,  then  t  =  1,  h  =  0,  en  =  0  for  n  i  0, 
f  =  1,  fn  =  0  for  nil,  and  (11)  reduces  to  (63),  of  course. 


TYPE  B,  (64) 

Let  v  =  +  •••  +  \>M  in  the  following.  The  asymptotic  decay 

of  (64)  then  agrees  with  (8).  Function  f(z)  in  (9)  is  given  by 


-  rT  (Y .  b ; :  .  b 

m=l  vr  m  n 


(75) 


while  product  function  F(z)  in  (68)  is 


f(z)  ■  t  TT  (i  -  v) 

m=l  ' 


(76) 


where 
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t  = 


(V  -  g°) 


H  <V-bma,” 


,  hm  = 


0  -  b 


m 


m  y  -  b  a 
m 


for  1  £  m  £  M  .  (77) 


In  order  to  develop  the  power  series  expansion  of  (76), 
observe  that 


M  <=° 

In  F ( z )  =  In  t  -  )  1  m  ln(l  -  hz)  =  )  e„  zn  , 

m=l  n=0 


(78) 


with 


e_  = 


In  t  ,  e„  =  -  u  , 
n  n  'n 


M 

u  =  )  vm  h*"}  for  nil 
rn  ‘ — Y  m  m 

m=l 


(79) 


Then,  by  use  of  (73)  -  (74),  we  find  coefficients 


n 


£o  =  £  '  £n  =  K  C  "k  £n-k  for  n  4  1  • 


(80) 


A  possible  initial  guess  for  p  follows  by  matching  the  n  =  0 
term  of  expansion  (11)  to  the  origin  behavior  of  the  given 
characteristic  function;  namely,  as  l  -»  0,  (11)  gives 


(1  -  Up)_v  ~  1  +  Hpv  , 


(81) 


while  (64)  yields 


M 


f (S)  ~  1  +  u 


m=l 


b  v 
m  m 


(82) 


Therefore,  we  could  choose  initially 


P  =  —  )  b  \)  ,  v  = 

\i  * - m  m  ' 


M 


M 


v  *— y  m  m 
m=l 


m=l 


m 


(83) 
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TYPE  C,  (65) 


Here,  the  asymptotic  decay  is  v  =  0  from  (65).  There  follows 


l  M  q  \ 

F(z)  =  f  (z)  =  exp  (a  -  z)  YU  i~-  ~h~  z 

v  m=l  m  ' 


(84) 


where 


= 


rm 


/  K  = 


P  -  c 


m 


*m  y  -  c  a  '  "  m  y  -  c  a 

1  m  1  m 


for  1  1  m  i  M  .  (85) 


Then, 


In  F(z)  =  (a  -  z)  ^  q  h£  zn  =  e  zn  ' 

m=l  m  n=0  m  n=0  n 


(86) 


where 


M 


=  a  >  '  q„  , 


M 


!n  “C  %  Cl  <ahm  -  D  tan!1’  <87> 
m=l 


m=l 


Then,  coefficients  ( f n )  follow  directly  from  (74). 

Sequence  |hm)  in  (85)  can  be  kept  relatively  small  by 

choosing  0  =  £  cm  /M.  Alternatively,  if  we  set  the  sum  of  jhm) 

to  zero,  we  find  6  =  E  c„  u  /£  u  ,  where  um  =  l/(y-acm).  On  the 

m  'm  ^m  rm  '  m' 

2  2  2 
other  hand,  minimizing  the  sum  of  (hi  gives  6  =  E  c„  u  /£  u  . 

3  ‘  m*  ^  m  rm  m 

Since  both  of  these  latter  two  approaches  give  a  0  value  which 
depends  on  y  and  a,  the  first  choice  seems  more  reasonable. 

A  smaller  { hm }  sequence  leads  to  smaller  jen)  and  (fn) 
sequences,  in  general,  thereby  promoting  faster  convergence  of 
expansion  (11). 
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TYPE  D,  (66) 


Again,  we  have  v  =  0.  Then,  we  find 


(  M  Qm 

F  ( z )  =  f(z)  =  exp  (y  -  pz)  YU  y~  h  z 

1  m=l  m 


(88) 


where 


m 


e  -  d, 


qm  y  -  d  a  ' 
m 


m 


hm  y  -  da 
m 


for  1  £  m  £  M 


(89) 


Then, 


M  00 

In  F  ( z )  =  (y  "  pz)  YU  <3, 


,  n  n  r — •  n 

,  h  z  =  )  e  z  , 
.  m  « — z  m  1 — *  n  ' 

m=l  n=0  n=0 


(90) 


where 


M 


M 


e  =  y  >  ;  q  ,  e^  =  >  ~  q_  h"  1  (yh  -  P)  for  nil.  (91) 

o  ‘ — y  m  n  c — y  m  m 

m=l  m=l 


Finally,  sequence  ( f n }  is  given  by  (74)  again. 

Sequence  (hm)  in  (89)  can  be  kept  relatively  small  by 
choosing 


P  = 


l 

M 


M 


m=l 


m 


(92) 


This  choice  is  independent  of  parameters  y  and  a. 
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APPLICATION  TO  FADING 


The  characteristic  function  of  the  output  of  a  multiple-pulse 
receiver  operating  in  a  noisy  medium  with  partially-correlated 
fading  was  derived  in  closed  form  in  [3;  (D-l)].  It  is 


Ul) 


(i  -  m)K(liM-1) 


x  exp 


M  K 

HZZZZ 

m=l  k=l 


emk 


(93) 


where  K  is  the  number  of  pulses,  M  is  the  number  of  fading 
components,  {Q^}  are  measures  of  random  signal-to-noise  ratios, 
and  { )  are  measures  of  deterministic  signal-to-noise  ratios. 

This  characteristic  function  in  (93)  is  recognized  to  be  the 
product  of  the  three  types  of  terms  A,  B,  C  in  (63),  (64),  (65), 
respectively.  This  means  that  the  logarithms  of  the  individual 
terms  will  add.  In  particular,  by  reference  to  (70),  it  means 
that  the  three  sets  of  (en)  coefficients,  of  each  type  of 
characteristic  function  component,  can  be  added  directly  together 
before  performing  the  exponentiation  indicated  by  (73)  and  (74). 
We  therefore  concentrate  on  obtaining  these  { en )  coefficients  for 
each  type  of  component  A,  B,  C  in  (93). 

With  regard  to  type  A  terms,  the  leading  term  in  (93)  leads 
us  to  identify  a  =  2,  v  =  -K(JjM-l)  in  (63).  Then,  by  reference 
to  (69)  and  (71),  we  find 
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=  (Y  -  ga) 

Iy  -  2aJ 


-K(ijM-l) 


h  = 


P  -  2 

Y  -  2a  ' 


■  t  .  en  -  -  £  (|  -  l)  (f^y"  <or,M 


(94) 


(95) 


The  fQ  coefficient  follows  from  (72)  as  f  =  t,  as  given  in  (94). 

As  for  type  B  terms,  the  second  term  in  (93)  requires  the 

identifications  M  ->  M, K;  m  ->  m,k;  b  ■>  20  ->  h  in  ( 64 ) . 

m  mk  m 

Then,  by  use  of  (83),  (77),  and  (79),  we  find  v  =  MK/2, 


t  = 


M  K 


(Y  -  Pa) 


MK/2 


(fr  TT(v  -  2aQmk)) 

k=l  mK  ' 


P  -  2Qmk 

h  -»  - 

m  y  ~  2aQ 


mk 


(96) 


e  =  In  t 
o 


,  M  K 

>  .  \ \ 

*n  ?n  *  *  ~  * 
n  m=l  k=l 


P  -  2Q 


mk 


Y  -  2aQ 


mk 


n 


for  nil 


(97) 


The  fQ  coefficient  is  given  by  (80)  as  fQ  =  t,  in  terms  of  (96). 

Finally,  for  type  C  terms,  the  last  term  in  (93)  yields  the 
identifications  M  -*  M,K;  m  ->  m,k;  ->  e  ■  ;  c  ->  2 Q  ■  in 
(65),  in  addition  to  v  =  0.  Then,  by  use  of  (85)  and  (87),  we 
find 


ink 

^m  y  -  2aQ 


h_  -» 


P  -  20 


mk 


mk 


m  y  -  2aQ 


mk 


(98) 


M  K  emR 

60  0  m=l  k=l  Y  ‘  2aQmk 


(99) 


n 


M  K  (g  -  2Q, 
=  (Pa  -  Y)  m  m 


mk 


n-1 


'mk 


m=l  k=l  (y  -  2aQmk) 


n+1 


for  nil 


(100) 
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The  f  coefficient  follows  from  (74)  in  the  form  fQ  =  exp(eQ), 
where  eQ  is  given  by  (99). 

When  we  combine  the  three  types  of  characteristic  function 
components  above,  the  three  individual  fQ  terms  multiply,  leading 
to  the  composite  value 


£  =  (Y  -  fSo)K 

O  ,  M  K 


(Y 


K(JjM-l) 


/  n  i\ 

TT  TT( 

vm=l  k=l 


Y  -  2aQ, 


mk 


) 


7 


exp 


M 


K 


a 


'mk 


m=t  k=T  ^  2aQ, 


mk 


(101) 


Similarly,  the  three  { en }  terms  given  by  (95),  (97),  and  (100) 
combine  (add)  to  yield  composite  value 


n  e  1  (P  -  2Qmk> 

n  2  m=l  k=l  (y  -  2aQmk) 


n-1 


n+1 


(«s  -  2Qmk)(Y  -  2aQmk )  + 


+  2n  *  V)  emR)  -  K(f  -  l)(e_5_yn  for  nil.  (102) 


The  composite  v  value  is  obtained  by  adding  the  three  values 
above,  to  get  v  =  K(1  -  M/2)  +  MK/2  +  0  =  K. 

The  coefficients  (f  1  can  be  found  from  (7^)  in  the  recursive 

form 

i  n 

f  =  -  y~Z  k  e.  f  .  for  nil.  (103) 

n  n  4— r  k  n-k  v  ' 


The  expansion  for  the  characteristic  function  is  again  given  by 
(11),  and  the  corresponding  probability  density  function  follows 
from  (19),  with  the  replacement  of  v  by  K  in  both  cases. 
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SPECIAL  CASE  OF  a  *  0 


In  this  special  case  of  a  =  0,  there  follows  f  =  1  from 
(101),  and 

Vn  n  en  =  i  ZZ  EZ  (ft  -  2Qmk)n"1  (P  -  2Qmk  -  2nemk)  + 

m= 1  R= 1 

+  K  (l  -  |)  (p  -  2 ) n  for  nil.  (104) 

The  characteristic  function  expansion  in  (11)  now  becomes 


f(S) 


1 

(1  -  ) K 


(105) 


Again,  we  could  write  fn(\>,0,O,y)  yn  =  ?n(v,P),  thereby 
emphasizing  the  two-parameter  dependence  in  this  special  case. 
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SPECIAL  CASE  OF  y  =  0 

Upon  use  of  y  =  0  in  (101)/  there  follows 


f 


o 


(106) 


Also ,  from  (102),  we  have 


n  i  M  K  (  R  »n-l 

an  n  e  =  ~  >  >  ( 1  -  — ] 

n  2  m=l  k=l  '  2Qmk^ 


_ 


2Q, 


mk 


1  -  n 


'mk 

^mk 


+  K(1  -  IK1  -  if  £or  n  *  1  • 


(107) 


The  characteristic  function  expansion  in  (11)  now  becomes 


HI) 


1 

(1  -  UP)K 


(108) 


Here,  also,  we  could  write  fn(v,0,a,O)  a11  =  fn(v,0),  explicitly 
indicating  the  two-parameter  dependence  in  this  special  case. 

As  a  check  on  the  above  results,  consider  the  situation  where 
there  is  no  signal  at  the  system  input  [3],  in  which  case  the 
parameters  in  (93)  become  e^  =  0  and  Qm^  =  1.  Then,  for  this 
noise-only  case,  (106)  yields  fQ  =  (0/2)  ,  and  (107)  reduces  to 

<xn  n  en  =  k|i  -  for  nil.  (109) 


The  corresponding  coefficients  { f n )  can  be  found  in  closed  form 
by  using  (73)  directly:  let  C  =  (1  -  0/2 )/a,  getting 
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00  00 

ezn  =  e  +  K  Y~2  ;  Cn  Zn  =  e  -  K  ln(l  -  Cz)  . 
n=0  n  °  n  ° 


(110) 


Then,  performing  the  exponentiation, 


-K 


(K) 


F(z)  =  exp(e  )  (1  -  Cz)  =  exp(e  )  } _ _ 

°  °  n=0  n' 


n  » n  n  . , . ,  . 

-  C  z  ,  ( 111 ) 


where  we  used  [8;  15.1.8].  There  follows  the  coefficients 


(K)  n  (K)„  1 1  «/o\n 

f  =  f  — ^  £,n  =  f  ^  1  -  P/2  for  n  *  0  . 

n  o  n!  ^  [2)  n!  [  a  ) 


(112) 


Then,  (62),  along  with  v  =  K,  yields  the  cumulative  distribution 
function 


for  threshold  v  >  0.  Since  this  result  must  be  independent  of  0, 
we  can  set  |3  =  2  and  get 

C ( v )  =  c(k,tJ)  for  v  >  0  .  (114) 

(This  checks  [3;  (188)  and  ( D— 17)]. )  The  independence  of  result 
(113)  on  £  has  been  numerically  verified. 
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APPLICATION  TO  NORMALIZATION 


The  general  results  derived  earlier  were  specialized  to  the 
case  where  parameter  y  =  0,  resulting  in  expression  (62)  for  the 
cumulative  distribution  function.  We  will  further  presume  that 
parameter  v  is  an  integer,  namely  K.  Then,  a  combination  of 
(62),  (13),  and  (B-14)  enables  us  to  write  the  corresponding 
exceedance  distribution  function  as 


E(v) 


C(v) 


)  ‘  f  a11  E  f K+n,^r|  for  v  >  0  . 

1^0  n  ~[ 


(115) 


Since  y  =  0,  we  let  f  (K,{S,a,0)  an  =  f  (K,0),  emphasizing  the 
two-parameter  dependence.  This  quantity  in  (115)  can  be 
interpreted  as  the  probability  that  a  random  variable,  denoted  g, 
exceeds  a  fixed  threshold  v: 


Pr( g  >  v)  =  E ( v )  » 


n=0 


f  E 
n  - 


K+n 


-I) 


for  v  >  0 


(116) 


Expansion  coefficients  ( f  n }  and  parameters  K  and  fi  are  arbitrary 
in  expansion  (116). 

Now,  let  us  suppose  that  random  variable  g  is  not  compared 
with  a  fixed  threshold,  but  rather  with  a  variable  threshold 
determined  from  some  finite  size  average  of  another  set  of 
independent  random  variables.  This  normalization  procedure  is 
opted  in  an  effort  to  realize  a  constant  false  alarm  receiver. 
In  particular,  let  random  variable  h  be  evaluated  from  the 
average  of  L  noise-cnly  squared  envelopes  outputs  of  disjoint 
narrowband  filters.  The  probability  density  function  of  h  is 
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then 


p ■„,»)  =  u--1-rxp^iu— ■' 

(2°)  (I.-l)! 


for  u  >  0  , 


(117) 


where  a  is  a  measure  of  the  noise  level  in  the  reference  bins. 

Random  variable  g  is  compared  with  a  scaled  version  of  random 
variable  h,  for  purposes  of  determining  whether  a  signal  is 
present  in  g.  This  leads  us  to  consider  the  exceedance 
(detection)  probability  of  the  normalizer  in  the  form 


Pr(g  >  v  h)  =  E ( v  h)  =  J  du  pQ(u)  E(v  u)  , 


(118) 


where  v  is  the  scale  factor.  Here,  we  used  (116)  and  (117),  and 
the  overbar  denotes  an  ensemble  average  over  random  variable  h. 
When  we  employ  (117),  (116),  and  (B-17),  probability  (118)  can  be 
developed  further  as  follows: 


Pr(g  >  v  h)  = 


f  du  .x£L-Wp_1  £  f  E [ K+n , 

A  (  2a )  L  (  L- 1 )  !  n=0  n  1 


xL  *  exp(-x/2)  ■ 


2  (  L- 1 )  !  n=0 


_ _  ?n  E  fK+n,^x]  = 

=  n  n  P 


f  dx 


xL  1  exp(-x/2 


2  ( L— 1 )  ! 


K-l+n 


_  £Xx  7^  L_ 

P  J  k!  1  Pxj 


K-l+n 


L-l+k 


=  YZ  £  ZZ  fr  ^  f  dx  -S -  exp  -x  i  +  SX  = 

K=o  n  k=0  k!  I  pj  J  2L  "l  U  |SJJ 
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which  can  be  easily  evaluated  by  recursions. 

Expansion  (119)  for  the  normalizer  detection  probability, 
Pr(g  >  v  h),  is  now  available  for  numerical  evaluation,  once  the 
expansion  coefficients  { f n }  in  (115)  have  been  determined.  No 
additional  integrals  need  to  be  evaluated;  rather,  only  the 
auxiliary  functions  in  (120)  must  be  evaluated.  This  procedure 
parallels  that  used  in  [4;  pages  13  -  16]. 
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EXTENSION  TO  ADDITIVE  GAUSSIAN  NOISE 


The  expansions  derived  above  are  most  applicable  to  the 
physical  situation  where  filter  outputs  have  been  envelope 
detected,  thereby  creating  positive  random  variables.  If  there 
is  also  present  additive  Gaussian  noise  at  the  system  output, 
modifications  to  the  expansions  must  be  expected  in  order  to 
accommodate  the  negative  values  that  are  then  possible  for  the 
random  variable  of  interest. 

Here,  we  will  derive  the  new  expansions  for  the  special  case 
where  parameter  y  is  taken  to  be  zero.  We  begin  by  repeating  the 
expansion  for  the  probability  density  function  in  (26): 


v— 1 


P^U) 


=  (U/P)v  exp(-u/g) 
fi  r(v) 


LI  £ 

n=0 


n 


ictuZil 

(v) 


n 


for  u  >  0 


(121) 


n 


It  is  presumed  that  this  expansion  has  already  been  accomplished 
for  random  variable  Xj.  A  second  random  variable,  which  is 

zero-mean  Gaussian,  is  added  to  x^.  The  probability  density 
function  of  is 

2 

pn ( u )  =  - exp f - for  all  u  .  (122) 

1  (2n)*o  K  2 aZ) 

The  actual  derivation  of  the  probability  density  function  of 
x  *  x^  +  *2  is  accomplished  in  appendix  C  by  convolving  results 
(121)  and  (122).  The  end  result  is 


p(u)  = 


rV  *  exp(-^u^/o^) 

PV  r( v) 


n=0 


n 


(gg/P) 

(v) 


n 


n 


[v-l+n  , 


u  o] 
o  "  pj 


(123) 
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for  all  u,  where 

00 

I/j(a)  5  J  dx  x^  (2n)-i*  exp|-  lx2  +  axj  = 

0 

=  a  ( a )  +  (p-1)  I^_2 (a)  for  p  >  1  .  (124) 

Starting  values  for  recurrence  (124)  when  p  is  integer,  that  is, 
when  v  in  (121)  is  integer,  are  given  in  appendix  C.  Thus,  the 
simple  recurrence  in  (124)  allows  for  rapid  evaluation  of  the 
probability  density  function  p(u)  in  (123),  once  coefficients 
j  f n }  in  (121)  have  been  determined. 

The  cumulative  distribution  function  corresponding  to  p(u)  is 
also  evaluated  in  appendix  C.  It  is  given  by 

C(v)  =  EZ  f„  a"  Jv_1+n£, f)  for  all  v  ,  (125) 

n=0  r 

where 


<ya,b) 


*  I  0*  *(a  -  f) 


(126) 


A  recurrence  for  is  derived  in  appendix  C  for  the  case  where  p 
is  integer;  in  this  case,  cumulative  distribution  function  C(v) 
in  (125)  can  be  rapidly  evaluated  once  the  coefficients  { f n }  are 
available.  Therefore,  the  addition  of  Gaussian  noise  is  easily 
accommodated,  in  terms  of  finding  the  density  and  distribution  of 
the  sum,  with  the  expansion  already  available  in  form  (121). 
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SUMMARY 

Expansions  of  probability  density  functions  and  their 
corresponding  cumulative  distribution  functions  have  been 
accomplished  in  a  form  that  can  be  interpreted  as  an  extension  of 
the  generalized  Laguerre  expansion.  This  form  does  not  require 
or  possess  any  orthogonality  properties.  Moments  or  cumulants  do 
not  have  to  be  extracted;  rather,  a  linear  fractional  transfor¬ 
mation  is  made  on  the  available  analytical  characteristic 
function,  followed  by  a  power  series  expansion.  Applications  to 
a  class  of  characteristic  functions  typically  encountered  in 
signal  processing  were  then  made  to  illustrate  the  procedure  by 
which  the  expansion  coefficients  can  be  obtained.  These  results 
are  exact,  not  approximations;  their  utility  depends  on  being 
able  to  determine  the  power  series  expansion  mentioned  above. 

These  results  have  already  been  used  to  advantage  in 
numerically  evaluating  the  performance  of  receivers  in  fading, 
both  with  and  without  normalization  [3,4].  In  some  cases, 
appropriate  choices  of  parameter  values  leads  to  series  of  all¬ 
positive  terms,  with  their  attendant  high  accuracy,  even  in  the 
range  of  very  small  probabilities  [3].  However,  it  must  be  noted 
that  these  results  are  very  example-dependent  and  may  not  extend 
simply,  in  general.  A  new  technique  and  expansion  will  usually 
be  required  for  each  different  problem. 

The  availability  to  choose  four  parameters,  v,p,a,Y*  allows 
for  better  fits  and  more  rapid  convergence  to  a  given 
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characteristic  function  or  probability  density  function.  This  is 
a  worthwhile  extension  of  the  Hermite  and  generalized  Laguerre 
expansions,  where  only  two  free  parameters  are  available.  Of 
course,  only  three  of  the  four  parameters  are  fundamental. 

Once  the  characteristic  function  expansion  for  a  fixed 
threshold  system  has  been  determined,  the  probability  density  and 
cumulative  distribution  functions  for  a  related  normalizer  follow 
immediately  without  the  need  for  extensive  new  derivations.  This 
observation  allows  for  a  ready  determination  of  the  cost  of 
employing  variable  detection  thresholds,  which  are  inherent  in  a 
constant  false  alarm  receiver. 
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APPENDIX  A.  FOURIER  TRANSFORM  OF  CUMULATIVE  DISTRIBUTION 

The  cumulative  distribution  function  C(v)  is  given  by 

v 

C(v)  =  J  du  p(u)  ,  (A-l) 

—  CO 

where  p(u)  is  the  probability  density  function,  and  f(£)  is 
the  characteristic  function  of  the  random  variable  under 
consideration.  The  two-sided  Fourier  transform  of  C(v)  is 

g( £ )  =  J  dv  exp(iSv)  C(v)  .  (A-2) 

This  latter  integral  converges  at  v  =  -»  if  <  a,  where  a  is 
a  positive  constant  (and  could  be  +°°,  as  for  a  positive  random 
variable).  The  integral  converges  at  v  =  +<»  if  ^  >  0.  Then, 
we  have 

v  +°° 

g(£)  =  J  dv  exp(i£v)  J  du  p(u)  =  J  du  p(u)  J  dv  exp(i£v)  = 

-OO  u 

=  J  du  p(u)  expt+H.),-  exp(Uu)  .  ^1  f(t)  t  (A-3> 

using  >  0.  That  is, 

9(0  =  t(D  for  0  <  li  <  a  ;  (A-4) 

this  relation  holds  within  a  strip  in  the  complex  £-plane. 

Let  contour  D  in  the  complex  S-plane  be  the  real  axis  with  a 
slight  upward  indentation  at  £,  =  0.  Then,  given  the 
characteristic  function  f(£),  we  can  return  to  the  cumulative 
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distribution  function  C(v)  by  means  of 


2^  |  exp(-ivS)  g(K)  =  C(v) 


(A-5) 


EXAMPLE  1:  TWO-SIDED  EXPONENTIAL 


p(u)  =  ^  exp(-|u|)  for  all  u 


(A-6) 


f(£)  =  - k  for  -1  <  $.  <  1  . 

i  +  r  1 


(A-7) 


Here,  a  =  1 .  Relation  ( A— 4 )  yields 


g(J.)  =  - = - for  0  <  l.  <  1  . 

uti  +  r ) 


(A-8) 


Then,  integral  (A-5)  yields,  with  contour  D  passing  below  l  =  i 
and  above  l  =  0,  the  correct  cumulative  distribution  function 


2^  |  exp( -ivl)  q(l)  = 


2  exp(v)  for  v  <  0 


1  -  exp(-v)  for  v  >  0 


=  C ( v ) .  ( A-9 ) 


EXAMPLE  2:  ONE-SIDED  EXPONENTIAL 


p(u)  =  exp(-u)  for  u  >  0  , 


(A- 10 ) 


C  ( v )  = 


1  -  exp(-v)  for  v  >  O' 


for  v  <  0, 


( A- 1 1 ) 
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f(Z)  =  (1  -  if;)-1  for  5.  >  -1  .  (A-12) 

Here,  a  =  +®.  From  (A-4), 

g(t>  -  for  o  <  ii  .  (A-13) 

Then,  (A-5)  yields  C(v)  as  given  by  (A-ll),  provided  that  contour 
D  passes  above  £  =  0. 

EXAMPLE  3 :  GAUSSIAN 


p(u)  =  (2n)-ii  exp(-u2/2)  for  all  u  ,  (A-14) 

v 

C(v)  =  4 ( v )  s  J  du  p(u)  for  all  v  ,  (A-15) 

-00 

f(Z)  =  exp( -S2/2 )  for  all  Z  .  (A-16) 

Here,  a  =  +°°.  Then,  from  (A-4), 

q(Z)  -  exp(-f;2/2)  for  0  <  .  ( A-17 ) 


With  contour  D  passing  above  Z  -  0,  we  find 

2 

2^  J  df;  exp( -iv£ )  gU)  =  J  dt  exp(-Hv)  expt~^  ^2) 
D  D 

=  -  ml exp(-iv*>  S*pt"p-2J  = 

ml  +  21  J d?;  ^Z^  exp<-*2/2>  = 
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=  \  +  ($(v)  -  =  *(v)  for  all  v  .  (A-18) 

Here,  we  used  [9;  3.896  4:  in  particular,  integrate  both  sides 
with  respect  to  b]. 
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APPENDIX  B.  INCOMPLETE  GAMMA  FUNCTION  RELATIONS 

Define  the  function 


A  V  1 

C(v,x)  =  J  dy  ^ ^  for  x  2  0  •  v  >  0  •  (B-l) 

0 

(This  is  P(v,x)  ^n  [8;  6.5.1].)  Then,  from  [8;  6.5.2  and  6.5.3], 


C(v,x) 


Y(v,x) 

r  ( v ) 


LL2Ll*1 
F(  v)  ' 


( B— 2 ) 


in  terms  of  incomplete  gamma  functions.  The  function  C(v,x)  is  a 
cumulative  distribution  function;  that  is, 

C(v,«)  =  1  .  ( B-3 ) 

It  is  also  useful  to  notice  that 


C(«,x)  =  0  ,  ( B-4 ) 

since  the  unit  area  of  the  integrand  of  (B-l)  moves  progressively 
farther  out  in  the  positive  direction  as  v  increases.  This  can 
be  seen  from  the  fact  that  the  density  (integrand)  in  (B-l)  has 
mean  v  and  standard  deviation  /v. 

Also,  define  (density)  function 

\i  —  1 

P(v,x)  =  * - r^P(~X)  for  x  2  0  ,  v  >  0  .  ( B-5 ) 

Then,  function  (B-l)  satisfies  recurrence 
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x 

C(v,x)  -  J  dy  p(v,y)  -  C(v-l,x)  -  p(v,x)  for  v  >  1  ,  (B-6) 

0 

while 

p(v,x)  -  p( v— 1  ,  x )  for  v  >  1  .  (B-7) 

Convenient  starting  values  for  (B-6)  and  (B-7),  when  v  is  integer 
or  half-integer,  are 

C(»j,x)  =  2  *  ( ( 2x  ) S)  -  1  ,  C(  1 ,  x )  =  1  -  exp(-x)  ,  (B-8) 

_  L 

p(h,x)  =  (nx)  7  exp(-x)  ,  p(l,x)  =  exp(-x)  .  (B-9) 

As  an  example,  with  v  =  n, 

p(n,x)  =  p(n-l,x)  —j-  for  n  =  2,3,4,...  (B-10) 

C(n,x)  =  C(n-l,x)  -  p(n,x)  for  n  =  2,3,4,...  ( B— 11) 

along  with  the  starting  values  for  p(l,x)  and  C(l,x)  given  in 
(B-9)  and  ‘B-8),  respectively 

As  a  second  example,  with  v  =  n  + 

p(n+*j,x)  =  p(n-*j,x)  for  n  =  1,2,3,...  (B-12) 

C(n+*5,x)  =  C(n-*j,x)  -  p(n+*j,x)  for  n  *  1,2,3,...  (B-13) 

along  with  the  starting  values  for  p(S,x)  and  C(*j,x)  given  in 
(B-9)  and  (B-8),  respectively. 


46 


TR  10371 


We  also  define  exceedance  distribution  function 


E(v,x)  =  1  -C{v,x)  =  J  dy  1 - rexp(-y) 

x 


=r(v,x)/r(v)  for  x  £  0  ,  v  >  0  . 

Then, 

E(v,x)  =  E(v-l,x)  +  p(v,x)  for  v  >  1  , 


with  recursion  (B-7)  for  p(\i,x),  and  starting  values 


E(Jj,x)  =  2  *|-(2x)J*j  ,  E ( 1 , x )  =  exp(-x)  . 

We  also  have  cloned  form 

n-1  k 

E(n,x)  =  exp(-x)  )  '  £7  for  n  £  1  . 

k=0  K" 


(B-14) 


(B-15 ) 


(B-16) 


( B—  1 7  ) 
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APPENDIX  C.  RECURRENCES  FOR  SUM  OF  RANDOM  VARIABLES 


The  probability  density  functions  of  random  variables  Xj  and 
X2  were  given  in  (121)  and  (122),  respectively.  We  will  now 
derive  the  probability  density  function  of  sum  x  =  +  x2,  ^y 

convolving  these  two  functions.  We  have 

P ( u )  =  |  dt  Pj(t)  p2(u-t)  = 


where  I  (a)  was  defined  in  (124).  If  v  is  integer,  the 
P 

recurrence  may  be  started  with 

IQ(a)  =  exp(a2/2 )  *(a)  ,  Ij(a)  =  a  IQ(a)  +  (2n)_1}  ,  (C-2) 

where  normalized  Gaussian  cumulative  distribution  function 


a 

*(a)  =  J  dx  <p(x)  ,  $(x)  =  (2n)~,?  exp(-x2/2)  .  (C-3) 


The  cumulative  distribution  function  corresponding  to  p(u)  is 

V  00 

C(v)  =  J  du  p ( u )  =  |  du  px(u)  C2(v-u)  = 

-00  0 


n=0 


'n  °"J 


du 


uv-l-Hn  exp^  _u/g ) 
pv+n  T(v+n) 
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,  n  _ 
n  a  v-l+n 


(5'f) 


for  all  v  , 


( C  —  4  ) 


where  J^(a,b)  was  defined  in  (126). 

A  recurrence  for  is  now  derived  for  the  case  where  p  is  an 
integer  m;  we  have,  for  b  >  0, 


Va'b> 


f  dx  *  e^i-x)  ,fa 

J  m ! 


-  §)  -  1  0*  - 


a-x/b 


exp( -x 
m! 


'  dt  <j>( t )  = 


b(a-t ) 


=  |  dt  4»  ( t )  J  dx  — 


exp( -x ’ 


=  |  dt  ♦(t) (l  - 


m  ,  k .  .  .  k 

— >  b  a-t 


ltl_  eb(t-a) 


m  uk  p  i 

-  #(a)  -  TZ  St  dt  *(t)  (a-t)K  exp(b( t-a)  )  = 
k=0  J 


m  ,  k 


)  -  II I  ur  f  dx  (2n)_i{  xk  expf-  i(a-x) 

.  n  n. .  J  z 


1  m  hk 

*(a)  -  exp(-a  /2)  YU  V7  1\c(a  "  b>  • 

k=0  K 


( C  —  5  ) 


There  follows  immediately  the  simple  recursion 


?  h 

Jm(a,b)  =  Jm_1(a,b)  -  exp(-a  /2)  im(a  -  b)  for  m  >  1  .  (C-6) 


This  can  be  started  with  the  value 


JQ(a,b)  =  $(a)  -  exp|^b2  -  abj  $(a  -  b) 


(C-7) 


60 


TR  10371 


Since  probability  density  function  exp(  -x)/r(/j+l )  in  (126) 
has  mean  p+1  and  standard  deviation  (/j+1)  ,  its  area  moves  to  the 
right  on  the  x-axis  as  p  increases.  Therefore,  the  $(a  -  x/b) 
function  in  (126)  forces  J^(a,b)  -»  0  as  p  ■*  ®.  This  observation 
implies  that  the  cumulative  distribution  function  expansion  in 
(125)  has  two  convergence  factors,  namely  f  an  and  Jv_i+n* 
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